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I.  INTRODUCTION 


Wraparound  Fins  (WAFs)  are  a  family  of  fins  that,  when  stowed,  conform  or  “wraparound” 
the  surface  of  a  cylindrical  body.  As  a  result  of  the  packaging  advantage  WAFs  have  over  planar 
fins,  WAFs  are  prevalent  on  tube-launched  missile  and  rocket  systems.  Several  fielded  missiles, 
rockets  and  munitions  utilize  WAFs  for  stability;  among  these  systems  are  Multiple  Launch 
Rocket  System  (MLRS),  The  Army  Combined  Arms  Weapon  System  (TACAWS),  Advanced 
Precision  Kill  Weapon  System  (APKWS),  Line-of-Sight  Antitank  (LOSAT),  Brilliant  Anti-Tank 
(BAT),  Compact  Kinetic  Energy  Missile  (CKEM),  Hydra-70,  and  variants  of  the  Zuni  rocket. 
Figure  1  shows  a  set  of  four  WAFs  both  stowed  around  the  body  of  a  rocket  and  deployed. 


Figure  1.  Packaging  Advantage  of  WAFs 


The  geometry  of  a  WAF  is  typically  determined  by  the  diameter  of  the  missile  and  the 
number  of  fins.  The  curved  span  of  the  WAF  is  typically  the  missile  circumference  divided  by 
the  number  of  fins,  and  the  angle  of  curvature  is  360  degrees  divided  by  the  number  of  fins.  A 
majority  of  the  systems  utilizing  WAFs  have  four  fins;  therefore,  a  WAF  with  90  degrees  of 
curvature  is  common.  However,  several  2.75-inch  rockets  are  equipped  with  three  WAFs  for 
stability.  BAT  employs  four  overlapping  WAFs  for  stability  with  a  curvature  angle  of 
180  degrees. 

WAFs,  however,  do  come  with  aerodynamic  peculiarities.  Systems  equipped  with  WAFs 
exhibit  significant  rolling  moments  at  zero  incidence.  The  “induced”  rolling  moment  is 
documented  as  a  function  of  Mach  number  and  Angle-of- Attack  (AoA). 

A.  Historical  Perspectives 

1.  United  States  Army 

A  series  of  tests  were  conducted  between  1971  and  1976  by  the  Aeroballistics 
Directorate  of  the  U.S.  Army  Missile  Research,  Development,  and  Engineering  Center 
(AMRDEC)  to  identify  alternative  stabilizing  devices  [6-10].  Among  these  devices  were  WAFs, 
ringtails,  and  flares.  Limited  data  were  collected  on  several  WAF  geometries  on  a  splitter-plate 
and  on  a  generic  4-inch  diameter  body  with  a  2-caliber  secant  ogive  nose  and  an  8-caliber 
cylindrical  after-body.  The  fins  tested  were  limited  to  90  degrees  of  curvature. 
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In  terms  of  stability,  the  U.S.  Army  concluded  that  WAFs  perform  similarly  to 
planar  fins  of  equivalent  projected  plan-form  shape.  It  was  also  noted  that  WAFs  produced  a 
substantial  amount  of  rolling  moment  which  varied  with  AoA  and  Mach  number.  These 
variations  in  rolling  moment  could  possibly  lead  to  significant  dynamic  problems  including 
Magnus  instability  and  roll  rate  variations  during  ballistic  flight  if  not  compensated  for  correctly. 
Furthermore,  the  rolling  moment  was  found  to  be  a  strong  function  of  Mach  number  as  the 
direction  of  the  rolling  moment  changed  near  Mach  1.0.  In  supersonic  flow,  the  fms  produced  an 
induced  normal  force  away  from  the  center  of  curvature  at  zero  incidence.  Conversely,  the  fins 
produced  an  induced  normal  force  toward  the  center  of  curvature  in  subsonic  flow  at  zero 
incidence. 


2.  United  States  Air  Force 

In  the  late  1980’s,  the  U.S.  Air  Force  (USAF)  began  investigating  the  cause  of  the 
low-incidence  rolling  moment  generated  by  their  tube-launched  missile  systems  equipped  with 
WAFs  [1,14,18].  The  USAF  used  several  techniques  to  investigate  the  flow  field  near  a  WAF 
including  free-flight  gun  tests,  wind-tunnel  tests  (with  and  without  the  aid  of  pressure-sensitive 
paint),  and  Computational  Fluid  Dynamics  (CFD).  The  USAF  also  investigated  several  methods 
of  reducing  the  magnitude  of  the  induced  roll  by  slotting  WAFs  and  altering  the  fin-body 
junction  angle.  A  majority  of  the  testing  was  performed  on  a  2.22  aspect  ratio  rectangular  fin 
with  a  thickness-to-chord  ratio  of  12.5  percent  and  a  45-degree  leading  edge  wedge  angle. 

Interest  was  focused  between  Mach  2.15  and  Mach  3.83. 

According  to  the  USAF  studies,  the  leading  edge  of  the  fin  causes  a  bow  shock 
that  interacts  with  the  convex  and  concave  sides  of  the  fin  much  differently.  On  the  concave  side 
of  the  fin,  the  shock  is  focused  near  the  center  of  curvature  causing  a  region  of  relatively  high 
pressure  which  diminishes  as  the  shock  becomes  more  acute  at  higher  Mach  numbers.  The 
convex  side  of  the  fin  shows  a  small  region  of  high  relative  pressure  near  the  body-fin  juncture 
that  intensifies  as  the  Mach  number  increases.  The  result  is  a  net  force  away  from  the  center  of 
curvature  which  decreases  with  Mach  number. 

B.  Current  Perspective 

The  U.S.  Army  Aviation  and  Missile  Command  tested  a  series  of  WAFs  on  a  splitter- 
plate  at  the  Lockheed  Martin  Missile  and  Fire  Control  High-Speed  Wind-Tunnel  (LMMFC 
HSWT)  in  Dallas,  Texas  in  January  of  2005,  with  the  goal  of  developing  a  design  methodology 
for  WAFs.  The  test  data  for  the  WAF  show  two  notable  features.  The  more  notable  feature  is  an 
induced  normal  force  on  the  WAF  at  zero  incidence  which  leads  to  an  induced  rolling  moment 
when  the  fins  are  used  on  a  missile  system.  The  second  difference  is  a  slight  increase  in  the 
normal  force  slope  with  respect  to  AoA  with  increasing  curvature.  Since  there  is  only  a  slight 
change  in  the  normal  force  slope,  it  appears  that  the  fin  curvature  effectively  generates  an 
induced  AoA  when  compared  to  a  planar  fin  of  the  same  projected  plan-form  shape. 
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II.  METHODOLOGY 


In  order  to  develop  a  design  methodology  for  WAFs,  the  effect  of  curvature  on  the  pressure 
loading  of  a  WAF  must  be  understood.  The  pressure-sensitive  paint  results,  presented  in 
Reference  14,  show  the  pressure  loading  of  a  WAF  at  zero  incidence  is  similar  to  the  pressure 
loading  of  a  planar  fin  at  an  AoA.  The  pressure  loading  has  distinct  divisions  that  appear  much 
like  Mach  lines.  The  interior  of  the  WAF  has  a  fairly  constant  pressure  and  the  tip  of  the  fin  has  a 
much  lower  pressure.  The  pressure  loading  is  similar  to  the  results  obtained  from  Evvard’s  theory 
for  a  planar  fin  at  a  non-zero  incidence.  Therefore,  it  is  reasonable  to  assume  that  the  pressure 
loading  of  a  WAF  can  be  estimated  with  Evvard’s  theory  with  the  addition  of  an  induced  AoA.  In 
addition  to  obtaining  the  normal  force  and  hinge  moment  of  the  fin,  the  geometry  of  the  WAF 
can  then  be  used  to  obtain  the  side  force  and  root  bending  moment  from  the  pressure  distribution. 

A.  Induced  AoA 

At  the  2005  LMMFC  HSWT,  fin-alone  data  was  gathered  via  a  splitter-plate  for  three 
different  aspect  ratio  rectangular  fins  with  various  curvature.  The  fins  were  attached  to  a  six- 
component  balance;  therefore,  a  complete  force-and-moment  data  set  was  gathered.  The  zero 
normal  force  AoA  of  each  tested  fin  was  derived  from  the  test  data,  and  a  correlation  dependent 
on  Mach  number,  aspect  ratio,  and  fin  curvature  was  formulated  for  the  induced  AoA.  The 
geometry  of  the  fins  tested  is  tabulated  in  Table  1 ,  and  a  photo  of  the  test  fins  can  be  seen  in 
Figure  2. 


Table  1.  Wind-Tunnel  Fin  Geometries 
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Figure  2.  Photo  of  Tested  WAFs 


B.  AoA  Dependence 

In  the  late  1940’s,  John  Evvard  [1 1,12]  and  others  [13,15,16]  solved  the  potential  flow 
equations  for  a  point-source  distribution  over  a  planar  fin  in  supersonic  flow.  In  order  to  utilize 
Evvard’s  solution,  the  fin  is  divided  into  regions  of  similar  disturbance  types  governed  by  the 
Mach  lines  emanating  from  leading  edge  discontinuities.  An  additional  region  can  form  on  swept 
fins  when  the  Mach  line  originating  from  the  root  leading  edge  discontinuity  is  reflected  by  the 
fin  tip  (Region  V  in  Fig.  3).  Each  region  consists  of  one  or  more  of  the  three  fundamental 
disturbance  types:  infinite  fin,  triangular  fin,  and  fin  tip.  The  potential  flow  solution  applicable  to 
each  region  is  used  to  determine  the  pressure  differential  of  the  upper  and  lower  surface  of  the  fin 
as  a  function  of  AoA. 

Since  the  regions  of  flow  are  defined  by  the  intersection  of  the  Mach  cones  and  the  fin 
surface,  curvature  can  have  a  significant  effect  on  the  zoning  of  the  fin  surface.  While  a  Mach 
cone  intersects  a  planar  fin  with  a  linear  Mach  line,  the  intersection  of  a  Mach  cone  and  a  WAF 
produces  a  curved  Mach  line.  As  the  curvature  increases,  the  area  of  the  fin  in  the  region  that 
creates  the  largest  pressure  differential,  Region  I,  also  increases.  The  result  is  an  increase  in  the 
normal  force  slope  with  respect  to  AoA  with  curvature.  Figure  3  illustrates  the  effect  of  curvature 
on  the  dividing  Mach  lines. 
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III.  THEORY 

The  theoretical  modifications  required  to  obtain  the  pressure  loading  on  a  WAF  surface 
begin  with  geometry.  In  order  to  apply  Evvard’s  theory,  the  fin  of  interest  must  be  divided  into 
incremental  surface  panels  with  a  control  point  in  the  center  of  each  panel.  The  curvature  angle 
and  projected  plan-form  fin  geometry  are  used  to  define  an  array  of  3-Dimensional  (3-D)  control 
points  and  the  local  surface  slopes  at  each  control  point.  The  fin  geometry  and  the  flow 
conditions  are  then  used  to  define  the  Mach  lines.  Once  the  control  points  are  zoned  based  on 
their  position  relative  to  the  Mach  lines,  Evvard’s  theory  is  used  to  determine  the  pressure 
differential  at  each  control  point.  Finally,  the  incremental  panel  area,  the  local  surface  slope,  and 
the  differential  pressure  coefficient  are  used  to  determine  the  normal  force,  hinge  moment,  side 
force,  and  root  bending  moment  coefficients  of  the  fin. 

A.  Curved  Fin  Geometry 

Defining  the  geometry  of  the  WAF  surface  is  the  basis  of  the  analysis.  The  fin  is 
divided  into  the  desired  number  of  span-wise  and  chord-wise  panels,  and  a  control  point  is 
positioned  in  the  center  of  each  control  panel.  With  the  chord-wise  (x)  and  span-wise  (y) 
coordinates  of  each  control  point  known,  the  magnitude  of  the  z-coordinate  is  determined  based 
on  the  curvature  of  the  fin.  Figure  4  shows  the  basic  nomenclature  that  will  be  used  to  describe 
the  geometry  of  a  WAF. 
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The  center  of  curvature  of  the  fin  is  defined  by: 

y,,  =  y™A.o  (l> 

*“=te) 

=  -4R'-yl  •  (3) 


Once  the  center  of  curvature  is  known,  the  z-component  of  the  fin  surface  can  be  obtained 
from  the  equation  of  a  circle  with  center  y0,  z0. 

z  =  Jr2  -{y-y0f  +  za  (4) 
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Furthermore,  the  local  surface  slope  of  each  control  point  will  be  used  to  obtain  the 
incremental  panel  area  on  which  the  pressure  differential  acts  to  produce  a  force  on  the  fin  in  the 
y-direction;  i.e.  side  force.  The  surface  slope  angle  at  each  control  point  is  defined  as: 


S  =  tan 


=  tan 


-i 


-(y-yj 

V*2 -{y-yof 


(5) 


B.  Dividing  Mach  Lines 

The  dividing  Mach  lines  of  a  WAF  are  derived  from  the  intersection  of  the  Mach  cone 
originating  at  the  fin  tips  and  the  fin  surface.  Figure  5  illustrates  the  intersection  of  the  two 
surfaces  showing  the  coordinates  that  are  referenced  in  Equations  (6)  through  (12). 


Figure  5.  Mach  Cone  -  WAF  Surface  Intersection 

For  a  planar  fin,  the  intersection  of  the  Mach  cone  emanating  from  the  fin  tip  and  the 
surface  is  simply  a  line  defined  by: 


tan  A  =  10'^  (6) 

where, 

A=sin"(X/)  •  (?) 
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As  seen  from  Figure  5,  the  line  describing  the  intersection  of  the  Mach  cone  and  a 
WAF  surface  can  be  redefined  to  include  the  z-component  as: 

The  Mach  cone  boundaries  and  their  reflection  lines  are  used  to  divide  the  fin  into  as 
many  as  five  regions  of  flow.  The  regions  shown  in  Figure  6  can  be  defined  as: 

Region  1 :  x  <  x/  and  x  <  X2 
Region  2:  x  >  x/  but  x  <  X2 
Region  3:  x  >  X2  but  x  <  x/ 

Region  4:x>  x /  and  x  >  xi  but  x  <  X3 
Region  5:  x  > 


y-axis 
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In  order  to  finalize  the  new  zoning  laws,  x/,  X2,  and  xj  must  be  defined  as  a  function  of 


y  and  z. 


(9) 

*1  =  P^y2  +  z2 

(10) 

x2  =  tanA  +  /?^/(l.0-y)2  +z 2 

(11) 

x,  =/3  +  /3yl(l  .0-y)2  +z2  . 

(12) 

The  Mach  cones  and  WAF  surface  intersections  are  represented  in  Figure  7  by  3-D 

surfaces. 


Figure  7.  3-D  Surface  Intersection 

In  order  to  validate  the  equations  used  to  zone  the  control  points  on  the  fin,  the  results 
at  Mach  1.6  for  a  rectangular  fin  with  a  chord  of  4.25  inches  and  a  span  of  3.0  inches  at  various 
angles  of  curvature  are  compared  to  the  3-D  Computer  Aided  Design  (CAD)  model.  Figures  8 
through  12  show  that  the  code  results  match  the  top-view  of  the  CAD  model  (Fig.  7)  for  various 
angles  of  curvature. 
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Figure  9.  Zoning  Verification  at  Mach  1.6  for  45.0  Degrees  of  Curvature 


Figure  10.  Zoning  Verification  at  Mach  1.6  for  90.0  Degrees  of  Curvature 
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Figure  11.  Zoning  Verification  at  Mach  1.6  for  135.0  Degrees  of  Curvature 


Figure  12.  Zoning  Verification  at  Mach  1.6  for  180.0  Degrees  of  Curvature 


C.  Pressure  Differential  in  Each  Region  of  Flow 

Now  that  the  fin  has  been  divided  into  zones  based  on  regions  of  influence,  the 
pressure  differential  between  the  upper  and  lower  surface  can  be  evaluated  based  on  the  types  of 
disturbances  that  affect  each  region  of  the  fin.  Since  the  potential  equation  for  a  fin  in  supersonic 
flow  is  described  by  an  ordinary  second-order  differential  equation,  the  laws  of  superposition 
apply.  Therefore,  the  pressure  differential  in  each  region  of  the  fin  is  a  summation  of  each 
upstream  disturbance  type.  Since  an  induced  AoA  method  is  being  utilized,  the  AoA  (a)  seen  in 
Equations  (14)  through  (25)  can  be  equated  to: 

^  =  ®  AERODYNAMIC  +  ®  INDUCED  •  (13) 


1  I 


1 .  Region  I 

Region  I  is  the  fundamental  portion  of  the  fin  which  lies  outside  both  Mach 
cones;  therefore,  control  points  within  Region  I  are  only  exposed  to  infinite  fin  (airfoil)  type 
disturbances.  From  linearized  supersonic  flow  theory,  the  pressure  coefficient  on  the  upper 
surface  of  a  flat  plate  is  given  by: 


and 


C 


p, lower 


2  a 

Vm2  -1 


C 


pjipper 


2  a 

Vm2-  i 


(14) 


(15) 


Differencing  the  lower  and  upper  pressure  coefficients  yield  a  differential  pressure  coefficient  of: 


A  C„ 


4  a 

Vm2-i 


(16) 


The  pressure  differential  in  Region  I  using  Evvard’s  theory  is  based  on  linearized  theory; 
however,  the  leading  edge  sweep  angle  is  included  such  that: 


4  a 


tan^A 


(17) 


2.  Region  II 


Region  II  is  located  within  the  interior  Mach  cone  that  is  produced  by  the  leading 
edge  discontinuity  at  the  root  of  a  swept  fin;  therefore,  it  is  referred  to  in  text  as  the  triangular  fin 
region.  Since  there  is  no  discontinuity  at  the  root  of  a  rectangular  fin  (A  =  0),  the  triangular  fin 
term  is  null,  and  Equation  (18)  reduces  to  Equation  (17).  Appropriately,  the  triangular  fin  effect 
increases  with  sweep  angle.  The  pressure  differential  in  Region  II  is  defined  as: 


where, 


AC 


Pji  ~ 


4a 

W  ~F2  —  tan2  A 


cos 


tan  A  +  pT_ 
TtanA  +  /? 


+  cos 


tan  A  -  f3T 
/J-TtanA 


(18) 


T  =  p^. 


X 


(19) 
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3.  Region  III 

Region  III  is  located  within  the  exterior  Mach  cone  that  is  produced  by  the 
leading  edge  fin  tip;  therefore,  it  is  referred  to  in  text  as  the  fin  tip  region.  Since  a  pressure 
differential  cannot  be  maintained  at  the  tip  of  a  fin,  the  potential  flow  equation  is  solved  with  a 
boundary  condition  imposed  such  that  the  pressure  differential  at  the  tip  of  the  fin  is  zero.  Region 
III  is  downstream  of  Region  I;  therefore,  the  tip  effect  is  an  addition  to  the  infinite  fin  solution. 
Since  the  tip  effect  uses  the  tip  of  the  fin  as  a  reference,  a  coordinate  system  is  defined  at  the 
leading  edge  fin  tip  such  that: 


xlip  =  x  -  tan  A 

(20) 

"5 

II 

1 

o 

(21) 

With  the  tip  coordinates  defined,  the  pressure  differential  coefficient  due  to  the  fin  tip 
disturbances  can  be  written  as: 


AC 


4a 


pytip 


- 


tan2  A 


COS 


-1 


xnp  +  y>,P(rP  +  tan  A)] 
x„p  -  y„r  tan  A 


(22) 


The  pressure  differential  in  Region  III  can  be  expressed  as: 


AC,.„,  =  AC,.,  +  AC, 


,///? 


(23) 


4.  Region  IV 

Region  IV  is  the  area  within  the  interior  and  exterior  Mach  cones;  therefore. 
Region  IV  is  affected  by  infinite  fin  disturbances,  triangular  fin  disturbances,  and  fin  tip 
disturbances.  Since  each  of  these  types  of  disturbances  have  been  defined,  the  pressure 
differential  in  Region  IV  is  simply: 


AC 


pJV 


A  Cpn  +  AC 


p,tip 


(24) 


5.  Region  V 


In  some  swept  fin  cases,  the  Mach  cone  originating  from  the  root  leading  edge 
discontinuity  intersects  the  fin  tip;  in  which  case,  an  addition  Mach  cone  is  created  with  an  origin 
at  the  fin  tip  intersection.  Thus,  the  fifth  fin  region  is  formed  within  Region  IV  designated  as 
Region  V.  Region  V  is  the  result  of  a  combination  of  Region  IV  disturbances  with  an  additional 
tip  effect  to  yield  a  pressure  differential  defined  as: 


^'p.V  ~ 


4a 


7r^j  01  —  tan2  A 


-1  xup  -  >v(tan  A  —  2/?)+ 2  tan  A 

COS  - —f - r - 

xnP  +  UV  +  2jtanA 


(25) 
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A  summary  of  the  pressure  coefficients  for  each  region  in  presented  in  Table  2. 


Table  2.  Pressure  Differential  Due  to  AoA 


Region 

Region 

Conditional 

Pressure  Coefficient  Differential 

I 

X  <  X/ 

and 

X  <  X2 

,n  4  a 

ACpj  r— - r— 

V/?"  -  tan"  A 

II 

X  >  Xi 

but 

X  <  X2 

4a 

a  r  -  - ™ 

M  01  -  tan2  A 

_)  tan  A  +  j8T  tan  A  -  BT 

cos  - - — hcos  - 

TtanA  +  p  y^-TtanA 

III 

X  >  X? 

but 

X  <  Xi 

AC,.,„=AC,,,+4C,# 

IV 

X  >  Xi 

and 

X  >  X2 

but 

X  <  Xi 

AC,.,,  =  AC„„  +AC„„ 

V 

X  >  Xi 

AZ-  _  . 4(X 

_ -i  X„P  ~  y,ip( tan  A  -  2/?)+  2  tan  A 

COS  - t - r - 

x„p  +  U,/,  +  2jtan  A 

D.  Empirically  Derived  Induced  AoA 

In  order  to  develop  an  empirical  expression  to  describe  the  induced  forces  and 
moments  generated  by  fin  curvature,  the  test  data  collected  at  the  January  2005  LMMFC  HSWT 
was  thoroughly  analyzed  to  find  a  correlation.  In  this  particular  test,  the  fins  were  mounted  on  a 
splitter-plate  to  minimize  the  appearance  of  shock  waves  upstream  of  the  fins.  Figure  13  shows 
one  of  the  WAFs  mounted  on  the  splitter-plate  along  with  the  test  sign  convention. 


14 


Figure  13.  WAF  on  Splitter-Plate  with  Sign  Convention 


In  order  to  obtain  an  empiral  relationship,  the  induced  AoA  of  the  three  different  aspect 
ratio  families  was  plotted  at  different  supersonic  Mach  numbers.  A  linear  curve-fit  was  used  to 
investigate  a  correlation  between  the  angle  of  curvature  and  the  induced  AoA.  Figures  14 
through  1 6  show  the  linear  relationship  of  the  three  aspect  ratio  fins. 
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INDUCED  ANGLE-OF-ATTACK  (DEGGREES)  INDUCED  ANGLE-OF-ATTACK  (DEGGREES) 


+  MACH  1 .5  ■  MACH  2.0  ♦  MACH  2.25  a  MACH  2.5  •  MACH  3.0 


Figure  14.  Curvature  Effects  for  AR  =  1.4118 


+  MACH  1 .5  ■  MACH  2.0  ♦  MACH  2.25  a  MACH  2.5  •  MACH  3.0 
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Figure  15.  Curvature  Effects  for  AR  =  1.8824 
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+  MACH1.5  ■  MACH  2.0  ♦  MACH  2.25  *  MACH  2.5  »  MACH  3.0  | 


Figure  16.  Curvature  Effects  for  AR  =  2.8333 


The  slopes  from  the  linear  fits  (Figs.  14  through  16)  were  then  plotted  with  respect  to 
aspect  ratio  at  each  Mach  number.  While  the  resulting  plots  varied  with  Mach  number,  the  aspect 
ratio  effects  were  matched  well  with  a  power  series  expression  for  each  Mach  number.  Figure  17 
shows  the  aspect  ratio  dependence  of  the  slope  at  Mach  2.25. 
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Up  to  this  point  in  the  analysis,  only  a  small  amount  of  the  collected  data  had  been 
analyzed  to  find  a  correlation.  Therefore,  the  next  step  was  to  include  all  the  data  in  the 
correlation.  The  results  from  the  test  data  were  tabulated  with  respect  to  curvature,  aspect  ratio, 
Mach  number,  and  induced  AoA.  In  order  to  capture  the  Mach  dependency,  a  genetic  algorithm 
was  used  to  find  three  different  power  series  relationships  that  could  be  applied  piecewise.  A 
dividing  Mach  number,  also  a  genetic  algorithm  variable,  would  be  used  to  capture  any  inflection 
in  the  data.  The  three  power  series  relationships  are  related  using  a  linear  interpolation  based  on 
Mach  number  about  the  dividing  Mach  number.  The  induced  angles  from  the  test  were  compared 
to  the  correlation,  and  the  genetic  algorithm  was  used  to  minimize  the  root  squared  sum  of  the 
differences  between  the  correlation  and  the  test  data.  A  second-order  polynomial  scheme  was 
also  investigated  with  less  success.  The  power  series  constants  and  dividing  Mach  number 
chosen  by  the  genetic  algorithm  and  the  associated  equations  are  presented  below.  In 
Equations  (26)  through  (30),  the  induced  AoA  ( conduced )  is  represented  in  degrees,  and  curvature 
angle  (6)  is  represented  in  radians. 

a  -  0.2425  •  6  •  ARum  (26) 

b  =  1.1583-  e-AR~2Am  (27) 

c  =  0.4346  •  9  •  AR~qa™  (28) 


For  Mach  >  2.2630 


a  INDUCED 


=  c  +  {b  -  c)- 


3.0  -Mach 
3.0-2.2630 


(29) 


For  Mach  <  2.2630 


a  INDUCED 


=  b  +  (a-b)- 


2.2630  -  Mach 
2.2630-1.5 


(30) 


As  seen  in  Figures  18  and  19,  the  correlation  worked  well.  The  induced  AoA  tends  to 
have  an  inflection  point  near  Mach  2.25.  The  dividing  Mach  number  helps  the  correlation  capture 
the  inflection  of  the  data.  The  overall  performance  of  the  correlation  is  seen  in  Figure  19.  While 
the  correlated  induced  AoA  is  not  quite  within  the  uncertainty  of  the  measured  induced  AoA 
(+-  0.1 -degree),  the  correlation  was  assumed  to  be  adequate  to  proceed  with  the  methodology. 
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Figure  18.  Correlated  Induced  AoA  Versus  Mach  Number 
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Figure  19.  Measured  Versus  Correlated  Induced  AoA 
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IV.  INTEGRATION  OF  THE  PRESSURE  DISTRIBUTION 


Using  the  induced  AoA  in  conjunction  with  the  modified  form  of  Evvard’s  theory  results  in 
a  method  of  predicting  the  induced  and  AoA  dependent  pressure  loading  on  a  WAF.  In  order  to 
obtain  the  force  and  moment  coefficients,  the  pressure  differentials  must  be  numerically 
integrated  over  the  WAF  surface. 

The  aerodynamic  loads  exerted  on  the  fin  are  a  product  of  the  differential  pressure 
coefficient  and  the  appropriate  fin  area.  The  normal  force  and  side  force  are  based  on  the  area  of 
the  fin  projected  onto  the  xy-plane  and  the  xz-plane,  respectively.  The  moments  are  a  summation 
of  the  product  of  the  incremental  forces  and  their  respective  distances  to  the  reference  location. 
The  root  bending  moment  is  referenced  about  the  root,  and  the  hinge  moment  is  referenced  about 
the  leading  edge.  The  forces  and  moments  can  be  expressed  as: 


Mv  v  =  Ar  •  Ay 


(31) 


(32) 


(33) 


(34) 


(35) 


S;>C/rAA.vo..-tan^. 


(36) 


(37) 
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V.  RESULTS 


The  results  of  applying  the  induced  AoA  to  Evvard’s  theory  modified  for  fin  curvature  show 
good  agreement  at  both  positive  and  negative  AoA.  The  results  are  presented  in  two  segments: 

(1)  selective  plots  comparing  theoretical  results  and  test  data  as  a  function  of  AoA  and,  (2) 
pressure  differential  contours  at  5  degrees  of  incidence  on  selective  fins.  The  normal  force,  side 
force,  root  bending  moment,  and  hinge  moment  will  be  compared  to  test  data  with  respect  to 
AoA.  Since  the  focus  is  on  the  pressure-driven  forces  and  moments,  the  axial  force  and 
associated  axial  moment  will  be  neglected  in  the  results. 

A.  Comparison  to  Test  Data 

Each  plot  contains  data  for  various  angles  of  curvature  from  both  the  methodology  and 
from  the  splitter-plate  test  data.  The  theoretical  solution  is  represented  as  a  line  (dotted,  dashed, 
and  solid).  The  test  results  are  represented  as  various  symbols  denoted  in  the  legend  on  each 
chart  along  with  the  word  “Test.”  The  test  data  represents  both  viscous  and  inviscid  effects,  while 
the  theoretical  solution  only  models  the  inviscid  effects.  One  source  of  discrepancy  between  the 
theoretical  solution  and  the  test  results  is  the  thickness  of  the  fins.  Theoretically,  the  fins  are 
infinitely  thin.  In  addition  to  a  finite  thickness,  the  test  fins  have  an  increased  thickness  near  the 
root  for  required  structural  properties  (Fig.  2).  The  test  data  uncertainty  as  quoted  by  LMMFC 
HSWT  can  be  seen  in  Table  3,  and  is  shown  on  select  coefficients  in  Figures  20  through  59.  In 
most  cases,  the  size  of  the  error  band  is  smaller  than  the  data  symbol. 


Table  3.  Coefficient  Uncertainty 


«CN 

«CY 

wCHM 

mCRBM 

±  0.0055 

±  0.0020 

±  0.0055 

±  0.0025 

1 .  Normal  Force 

Figures  20  through  29  show  the  comparison  of  the  theoretical  normal  force  and 
the  test  results.  Overall,  the  addition  of  the  curvature  term  in  the  zoning  laws  allows  Evvard’s 
theory  to  model  the  changes  in  the  normal  force  slope  extremely  well  (within  the  uncertainty  of 
the  test  data).  The  normal  force  slope  of  the  moderate  aspect  ratio  (AR  =  1.41 18  and  AR  = 
1.8824)  fins  are  modeled  slightly  better  than  the  extreme  aspect  ratio  fins.  For  the  higher  aspect 
ratio  cases  (AR  =  2.8333)  (Figs.  26  through  28)  the  normal  force  slope  appears  to  be  slightly 
different  at  positive  versus  negative  AoA.  Evvard’s  theory  works  well  at  Mach  1.5,  but  the 
accuracy  seems  to  improve  with  increasing  Mach  number.  As  for  the  induced  normal  force,  the 
empirical  fit  for  the  induced  AoA  allows  the  normal  force  to  match  the  test  data  at  low  AoA.  The 
largest  discrepancy  in  the  low  incidence  normal  force  modeling  is  seen  in  Figure  20  for  a  fin 
curvature  of  180  degrees.  Even  though  leading  edge  sweep  angle  was  not  accounted  for  in  the 
empirical  fit,  Figure  29  shows  good  agreement  between  theory  and  test  data  for  a  leading  edge 
sweep  angle  of  35  degrees. 
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Figure  20.  Normal  Force  Comparison  for  AR  =  1.4118  at  Mach  1.5 


Figure  21.  Normal  Force  Comparison  for  AR  =  1.4118  at  Mach  2.25 
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Figure  22.  Normal  Force  Comparison  for  AR  =  1.4118  at  Mach  3.0 
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Figure  24.  Normal  Force  Comparison  for  AR  =  1.8824  at  Mach  2.25 


Figure  25.  Normal  Force  Comparison  for  AR  =  1.8824  at  Mach  3.0 
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Figure  26.  Normal  Force  Comparison  for  AR  =  2.8333  at  Mach  1.5 
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Figure  28.  Normal  Force  Comparison  for  AR  =  2.8333  at  Mach  3. 0 
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2.  Side  Force 


The  side  force  of  the  fin,  as  seen  in  Figure  13,  is  the  force  that  is  directed  from  the 
fin  tip  to  the  root.  The  pressure  contour  plots  (Figs.  60  through  71)  can  be  used  to  visualize  the 
manner  in  which  the  pressure  exerts  a  side  force  on  the  fin.  A  curved  fin  has  surface  area  in  the 
xz-plane;  therefore,  the  pressure  differential  generates  a  side  force  on  a  WAF. 

The  theoretical  solution  for  all  fins  and  conditions  shows  a  pushing  (directed  from 
tip  to  root)  side  force  at  positive  AoA  and  a  pulling  (directed  from  root  to  tip)  side  force  at 
negative  AoA.  The  test  data  shows  that  a  pulling  force  is  dominant  even  at  positive  AoA.  As 
seen  in  Figure  2,  the  tested  fins  were  thickened  near  the  root  for  structural  rigidity,  which  alters 
the  area  of  the  fin  in  the  xz-plane.  In  addition  to  the  thickened  root,  the  splitter-plate  also  is 
suspected  to  influence  the  side  force.  For  finite  thickness  fin  in  supersonic  flow,  the  leading  edge 
of  the  fins  will  produce  a  shock  wave.  On  the  convex  side  of  the  tested  WAFs  near  root,  the 
shock  wave  is  most  likely  reflected  by  the  splitter-plate  forming  a  region  of  high  pressure.  A 
region  of  high  pressure  in  this  location  would  generate  a  pulling  side  force.  Again,  the  theory 
does  not  account  for  fin  thickness.  However,  the  side  force  order  of  magnitude  and  data  trend  is 
modeled  well.  Figures  30  through  39  show  the  side  force  coefficient  comparison  of  the  theory  to 
the  test  data.  The  theory  typically  over-predicts  the  magnitude  of  side  force  at  positive  AoA,  but 
matches  the  magnitude  well  at  negative  AoA.  The  computed  side  force  of  the  swept  fin  matches 
test  data  remarkably  better  than  the  rectangular  fins  (Fig.  39). 


Figure  30.  Side  Force  Comparison  for  AR  =  1.4118  at  Mach  1.5 
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Figure  31.  Side  Force  Comparison  for  AR  ~  1.4118  at  Mach  2.25 


Figure  32.  Side  Force  Comparison  for  AR  =  1.4118  at  Mach  3.0 
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Figure  33.  Side  Force  Comparison  for  AR  =  1.8824  at  Mach  1.5 


Figure  34.  Side  Force  Comparison  for  AR  =  1.8824  at  Mach  2.25 
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Figure  35.  Side  Force  Comparison  for  AR  =  1.8824  at  Mach  3.0 


Figure  36.  Side  Force  Comparison  for  AR  =  2.8333  at  Mach  1.5 
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Figure  38.  Side  Force  Comparison  for  AR  =  2.8333  at  Mach  3.0 
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Figure  39.  Side  Force  for  AR  =  1.8824;  A  =  35°  at  Mach  1.5 
3.  Root  Bending  Moment 

The  span-wise  CP  of  the  fins  is  modeled  well  with  the  theory  as  reflected  in  the 
root  bending  moment  comparisons  seen  in  Figures  40  through  49.  The  root  bending  moment  is  a 
combination  of  the  normal  and  side  forces  with  their  respective  moment  arms  (Equation  (37)). 
While  the  root  bending  moment  is  modeled  accurately  at  low  incidence,  the  theory  tends  to 
under-predict  the  root  bending  moment  at  higher  AoA.  The  root  thickness  is  most  likely  moving 
the  span-wise  CP  of  the  tested  fins  towards  the  fin  tip,  thus  increasing  the  root  bending  moment. 
At  negative  AoA,  the  root  bending  moment  is  less  sensitive  to  fin  curvature. 
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Figure  40.  Root  Bending  Moment  Comparison  for  AR  =  1.41 18  at  Mach  1.5 
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Figure  41.  Root  Bending  Moment  Comparison  for  AR  =  1.4118  at  Mach  2.25 
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Figure  42.  Root  Bending  Moment  Comparison  for  AR  =  1.41 18  at  Mach  3.0 


Figure  43.  Root  Bending  Moment  Comparison  for  AR  =  1.8824  at  Mach  1.5 
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Figure  44.  Root  Bending  Moment  Comparison  for  AR  =  1.8824  at  Mach  2.25 


Figure  45.  Root  Bending  Moment  Comparison  for  AR  =  1.8824  at  Mach  3.0 
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Figure  46.  Root  Bending  Moment  Comparison  for  AR  =  2.8333  at  Mach  1.5 


Figure  47.  Root  Bending  Moment  Comparison  for  AR  =  2.8333  at  Mach  2.25 
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Figure  48.  Root  Bending  Moment  Comparison  for  AR  =  2.8333  at  Mach  3.0 


Figure  49.  Root  Bending  Moment  for  AR  =  1.8824;  A  =  35°  at  Mach  1.5 
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4.  Hinge  Moment 


The  hinge  moment  data  are  presented  in  Figures  50  through  59  about  the  mid-root 
chord  since  fins  are  typically  hinged  about  their  theoretical  Center-of-Pressure  (CP).  As  a  result, 
the  hinge  moments  are  very  small,  and  the  differences  between  theoretical  and  test  results  are 
visually  amplified.  Presenting  the  data  about  the  mid-root  chord  gives  a  better  indication  of  the 
chord-wise  CP  accuracy.  The  test  results  indicate  that  the  CP  is  more  forward  than  the  theoretical 
results.  Since  the  test  articles  were  not  infinitely  thin,  a  shock  wave  from  the  leading  edge  is 
suspected  to  decrease  the  Mach  number  over  the  fin,  which  would  shift  the  CP  forward.  As  seen 
in  Table  4,  the  chord-wise  CP  of  the  fin  is  modeled  within  5-  to  10-percent  relative  error  of  the 
test  results.  As  a  result  of  matching  both  the  normal  force  coefficient  and  the  chord-wise  CP 
reasonably  well,  the  hinge  moment  comparison  yields  reasonable  results. 


Table  4.  Chord-Wise  CP  Non-Dimensionalized  by  Lref 


Wind-Tunnel 

Theory 

Percent 

Wind-Tunnel 

Theory 

Percent 

Wind-Tunnel 

Theoiy 

Percent 

AR 

e 

Mach  1.5 

Mach  1.5 

Difference 

Mach  2.25 

Mach  2.25 

Difference 

Mach  3.0 

Mach  3.0 

Difference 

1.4118 

0 

0.3833 

0.4265 

10.12% 

0.4398 

0.4663 

5.69% 

0.4383 

0.4786 

8.41% 

1.4118 

45 

0.3845 

0.4264 

9.81% 

0.4413 

0.4663 

5.37% 

0.4786 

** 

1.4118 

90 

0.3951 

0.4257 

7.18% 

0.4383 

0.4664 

6.04% 

0.4425 

0.4791 

7.66% 

1.4118 

135 

0.3892 

0.4241 

8.21% 

0.4428 

0.4674 

5.25% 

0.4638 

0.4812 

3.60% 

1.4118 

180 

0.3914 

0.4215 

7.13% 

0.4487 

0.4724 

5.02% 

0.4655 

0.4869 

4.40% 

1.8824 

0 

0.4081 

0.4487 

9.04% 

0.4408 

0.4770 

7.590/0 

** 

0.4854 

** 

1.8824 

90 

0.4122 

0.4482 

8.03% 

0.4589 

0.4775 

3.90% 

0.4496 

0.4862 

7.53% 

2.8333 

0 

0.4207 

0.4706 

10.59% 

0.4391 

0.4867 

9.790/0 

** 

0.4918 

** 

2.8333 

90 

0.4158 

0.4708 

11.68% 

0.4409 

0.4875 

9.56% 

0.4200 

0.4926 

14.750/o 

2.8333 

135 

0.4309 

0.4722 

8.76% 

0.4610 

0.4897 

5.880/0 

0.4675 

0.4946 

5.49% 
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Figure  50.  Hinge  Moment  About  Cr/2.0  Comparison  for  AR  =  1.4118  at  Mach  1.5 
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Figure  51.  Hinge  Moment  About  Cr/2.0  Comparison  for  AR  =  1.4118  at  Mach  2.25 
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Figure  52.  Hinge  Moment  About  Cr/2.  0  Comparison  for  AR  =  1.4118  at  Mach  3. 0 


Figure  53.  Hinge  Moment  About  Cr/2.0  Comparison  for  AR  =  1.8824  at  Mach  1.5 


40 


0.04 


S 

x 


U 


0 


-0.04 


-10  0  10 

Angle-of- Attack  (degrees) 

Figure  54.  Hinge  Moment  About  Cr/2.0  Comparison  for  AR  =  1.8824  at  Mach  2.25 
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Figure  55.  Hinge  Moment  About  Cr/2.  0  Comparison  for  AR  =  1.8824  at  Mach  3. 0 
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Figure  56.  Hinge  Moment  About  Cr/2.0  Comparison  for  AR  =  2.8333  at  Mach  1.5 
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Figure  5  7.  Hinge  Moment  About  Cr/2.  0  Comparison  for  AR  =  2. 8333  at  Mach  2. 25 
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Figure  58.  Hinge  Moment  About  CRJ2.0  Comparison  for  AR  =  2.8333  at  Mach  3.0 


Figure  59.  Hinge  Moment  About  Cr/2.0  for  AR  =  1.8824;  A  =  35°  at  Mach  1.5 
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B.  Pressure  Contour  Plots 


Pressure  contour  plots  are  presented  in  Figures  60  through  71  for  5  degrees  AoA  to 
show  the  effects  of  fin  curvature  and  Mach  number.  The  pressure  contour  plots  can  give  another 
dimension  to  the  aerodynamic  assessment  of  WAFs.  The  CPs  can  be  visualized  and  areas  of 
high-pressure  differential  can  be  identified  for  possible  fin  redesign.  The  figures  are  arranged  to 
show  the  effects  of  Mach  number  on  each  row  and  the  effects  of  fin  curvature  on  each  column. 

The  CP  moves  aft  and  toward  the  fin  tip  as  the  Mach  number  increases.  The  area  of 
higher  pressure  differential  (Region  I)  enlarges  with  increasing  fin  curvature.  Potential 
performance  enhancements  could  be  generated  based  on  analysis  of  the  pressure  contour  plots 
such  as  clipping  the  trailing  edge  fin  tip  to  decrease  the  fin  surface  area,  while  retaining  the  fin 
region  which  produces  a  majority  of  the  stabilizing  force. 

Fin  curvature  also  increases  the  area  of  the  fin  in  the  xz-plane,  which  creates  a  higher 
side  force.  Since  fins  are  typically  used  in  sets,  the  net  fin  side  forces  are  cancelled;  however,  the 
side  force  can  generate  a  roiling  moment  when  the  CP  is  located  at  a  non-zero  z-coordinate. 
While  the  induced  fin  normal  force  generates  a  majority  of  the  induced  rolling  moment  on  a 
missile  equipped  with  WAFs,  the  fin  side  force  contributions  must  be  known  to  accurately  model 
the  overall  missile  rolling  moment. 
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Figure  60.  Pressure  Contour  for  0  =  0.0; 
AR  =  1.4118  at  Mach  1.5 
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Figure  61.  Pressure  Contour  for 
0  =  0.0;  AR  =  1.4118  at  Mach  3.0 


Figure  62.  Pressure  Contour  for  6  =  90.0; 
AR  =  1.4118  at  Mach  1.5 


Figure  63.  Pressure  Contour  for 
0  =  90.0;  AR  =  1.4118  at  Mach  3.0 


45 


Figure  64.  Pressure  Contour  for  6  =  180.0; 
AR  =  1.41 18  at  Mach  1.5 


Figure  65.  Pressure  Contour  for 
0  =  180.0;  AR  =  1.4118  at  Mach  3.0 


Figure  66.  Pressure  Contour  for  A  =  35. 0; 
0  =  0.0;  AR  =  1.8824  at  Mach  1.5 


Figure  67.  Pressure  Contour  for 
A  =  35.0;  0  =  0.0;  AR  =  1.8824  at 
Mach  3. 0 
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Figure  68.  Pressure  Contour  for  A  =  35. 0; 
0  =  90.0;  AR  =  1.8824  at  Mach  1.5 


Figure  69.  Pressure  Contour  for 
A=  35.0;  0  =  90.0;  AR  =  1.8824  at 
Mach  3. 0 
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Figure  70.  Pressure  Contour  for  A  =  35. 0; 
0  =  180.0;  AR  =  1.8824  at  Mach  1.5 
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Figure  71.  Pressure  Contour  for 
A  =  35.0;  0  =  180.0;  AR  =  1.8824  at 
Mach  3. 0 


VI.  LIMITATIONS 


While  several  fin  parameters  and  flow  conditions  were  used  to  present  a  generalized  method 
of  estimating  the  pressure-driven  forces  and  moments  of  a  WAF,  the  methodology  has  Mach 
number,  aspect  ratio,  leading  edge  sweep  angle,  and  thickness-to-chord  ratio  limits  based  on  the 
tested  fin  parameters.  The  empirical  fit  for  induced  AoA  is  based  on  data  collected  on  rectangular 
fins  from  Mach  1.5  to  Mach  3.0,  with  aspect  ratios  ranging  from  1.4 11 8  to  2.8333.  The 
thickness-to-chord  ratio  for  all  the  fins  varied  linearly  from  3  percent  at  the  root  to 
1 .5  percent  at  the  tip  (with  the  exception  of  the  thickened  root).  The  method  compares  well  to 
wind-tunnel  results  for  these  fins;  however,  no  proof  exists  that  this  correlation  applies  to  fins 
with  features  outside  the  test  envelope.  The  methodology  to  compute  the  pressure  differential  is 
derived  from  potential  flow  theory  for  an  infinitely  thin  fin;  therefore,  the  results  from  the 
method  will  represent  a  thin  fin  with  no  boundaiy  layer  or  shock  waves  that  will  never  stall.  With 
these  limitations  stated,  the  method  will  provide  a  reasonable  design  envelop  for  typical  WAF 
applications  in  supersonic  flow. 

VII.  CONCLUSION 

A  new  method  has  been  developed  to  obtain  the  pressure  loading  of  curved  or  WAFs  with 
modifications  to  existing  supersonic  fin  theory  and  the  aid  of  recent  wind-tunnel  test  data.  The 
theoretical  results  show  agreement  to  wind-tunnel  test  data  for  the  low-incidence  forces  and 
moments  that  have  puzzled  aerodynamicist  for  years.  The  method  provides  an  expedient  analysis 
tool  that  has  been  developed  to  cover  a  broad  range  of  fin  parameters  that  can  be  used  for  roll 
tailoring  or  roll  minimization  on  missile  systems  requiring  the  use  of  WAFs. 
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LIST  OF  ACRONYMS  AND  SYMBOLS 


AoA 

Angle-of-Attack 

AMRDEC 

Army  Missile  Research,  Development  and  Engineering  Center 

APKWS 

Advanced  Precision  Kill  Weapon  System 

BAT 

Brilliant  Anti-armor  Technology 

CAD 

Computer-Aided  Drafting 

CFD 

Computational  Fluid  Dynamics 

CP 

Center-of-Pressure 

CKEM 

Compact  Kinetic  Energy  Missile 

HSWT 

High  Speed  Wind-Tunnel 

LMMFC 

Lockheed  Martin  Missile  and  Fire  Control 

LOSAT 

Line-of-Sight  Antitank 

MLRS 

Multiple  Launch  Rocket  System 

TACAWS 

The  Army  Combined  Arms  Weapons  System 

USAF 

United  States  Air  Force 

WAF(s) 

Wraparound  Fin(s) 

Cn 

Normal  Force  Coefficient 

Crbm 

Root  Bending  Moment  Coefficient 

Cy 

Side  Force  Coefficient 

Chm 

Hinge  Moment  Coefficient 

&Ax.y 

Incremental  Fin  Panel  Surface  Area  Projected  onto  xy-Plane 

Ax 

Incremental  Chord-wise  Length 

Ay 

Incremental  Span-wise  Length 

Az 

Incremental  Curve-wise  Length 

b 

Wingspan 

b/2 

Fin  semi-span 

Lref 

Mean  Chord  Length 

Sref 

Fin  Plan-form  Area,  Lref  (b/2) 

AR 

Fin  Aspect  Ratio,  b2/(2.0  Sref) 

a 

Angle-of-Attack 

CX INDUCED 

Induced  Angle-of-Attack 

Ot, AERODYNAMIC 

Angle  Between  the  Free-stream  Mach  number  and  the  Fin  Chord 

M,  Moo 

Free-stream  Mach  number 

b 

Compressibility  Factor,  VM2  -1.0 

5 

Fin  Curvature  slope  angle 

A 

Leading  Edge  Sweepback  Angle 

6 

Fin  Curvature 
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